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Abstract: We calculate the collisional energy loss of a muon of high energy E in 
a hot QED plasma beyond logarithmic accuracy, i.e., we determine the constant 
terms of order O {1) in —dE/dx oc Ini? + 0{1). Considering first the t-channel 
contribution to —dE/dx, we show that the terms ~ O {1) are sensitive to the full 
kinematic region for the momentum exchange q in elastic scattering, including large 
values g ~ C [E) . We thus redress a previous calculation by Braaten and Thoma, 
which assumed q <^ E and could not find the correct constant (in the large E limit). 
The relevance of 'very hard' momentum transfers then requires, for consistency, that 
s and u-channel contributions from Compton scattering must be included, bringing 
a second modification to the Braaten-Thoma result. Most importantly, Compton 
scattering yields an additional large logarithm in —dE/dx. Our results might have 
implications in the QCD case of parton collisional energy loss in a quark gluon 
plasma. 
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1. Introduction 

Jet quenching, as anticipated by Bjorken 25 years ago [1], is a prominent signature 
of the intriguing state of matter created at the Relativistic Heavy Ion Collider. The 
suppression of light hadron spectra at large transverse momentum p± [2, 3] can be 
explained - at least qualitatively - by attributing the hadron attenuation to the 
radiative energy loss of the parent parton (light quark or gluon) induced by its 
rescatterings in the hot or dense medium. On the other hand, recent experimental 
data on heavy flavour quenching [4, 5] , measured indirectly via the p_L-spectra of 
electrons from D and B meson decays, suggest that the radiative energy loss of heavy 
quarks might be insufficient to explain the observed attenuation^. This renewed the 
interest in the collisional part —AEcoii of the parton energy loss [7], which in the 

^It is stressed in Ref. [6] that this statement might be somewhat premature, since the theoretical 
calculation of heavy quark production suffers from large uncertainties already in proton-proton 
collisions, and also because the contributions to the electron spectra from charm and beauty are 
not separated experimentally. 
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case of a heavy quark might not be neghgible - contrary to what has often been 
assumed. A basic quantity required to estimate colhsional quenching is the rate of 
energy loss per unit distance, —dE/dx, of a parton produced in the remote past 
and travelhng in a large size medium, as studied in Refs. [1,8,9]. For heavy ion 
collisions, where a parton initially produced in a hard subprocess crosses a medium 
of finite size L, we expect deviations from the linear law —AEcoii{L) = {—dE/dx) ■ L 
[10-14]. However, the knowledge of —dE/dx is a prerequisite before attempting any 
evaluation of —AEcoii- 

So far, the most detailed calculation of —dE/dx for a heavy quark in the quark 
gluon plasma is done by Braaten and Thoma [9], and is based on their previous 
evaluation of muon colhsional energy loss in QED [15]. Here we will reconsider 
the latter calculation, which appears to suffer from an incorrect assumption on the 
magnitude of the momentum exchange in elastic scattering. As in Ref. [15] (referred 
to as BT in the following) we study the propagation of a muon of mass M and 
momentum P = {E,p) in an e^7 plasma at a temperature T <^ M, but large 
enough to neglect the electron mass. The muon can be considered as a test particle, 
losing (or gaining) energy due to elastic (2 —* 2) scattering off thermal particles of 
momentum K. The latter can be electrons and positrons (t-channel scattering, see 
Fig. [|a), or thermal photons (Compton scattering, see Fig. The Mandelstam 
invariants of the elastic processes are defined as 

s={P + K)\ t={P-P')\ u = {P-K'f . (1.1) 

In the high energy limit^, the BT calculation incorrectly assumes that the (spacelike) 
momentum exchange P — P' = Q = [u , q) is always small compared to the incoming 
muon energy, namely \uj\ < q = \q\ <^ E. We mention that this problem has also 
been noted recently in the QCD context of Refs. [16, 17]^. 

As a first consequence of the BT assumption q <^ E, the Compton contribution 
—dEj/dx was neglected in Ref. [15]. Second, in the f — > 1 limit the result for the 
energy loss -dE^/dx due to scattering off electrons (and positrons) found in [15], 



dE, 



dx 



487r 



In ^ + 2.031 



;i.2) 



has an incorrect constant next to the leading logarithm. For a heuristic argument 
showing the importance of large momentum exchange q, we note that the logarithmic 



2 When E becomes large, s = IVP + 2PK ~ 2PK ~ O (ET), since K < O (T) is constrained 
by a thermal distribution. The high energy limit s ^ is thus equivalent to ^ M"^ /T. For 
convenience it will often be referred to as the v = pj E —f \ limit. 

^Those references indeed mention the need for a careful treatment of the kinematics, including the 
region of large transfers g ~ O {E). However, those studies seem to focus on t-channel scattering. 
In the present paper we also stress (in the case of QED), that including Compton scattering is 
required for consistency. 
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energy dependence in ( p. .21) arises from an integral J^™'"' dq/q, with g^ax — -E at high 
energy. While the region q ~ gmax does not contribute to the leading logarithm, 
it does contribute to the constant next to it {e.g., the interval [^max/S, Q'max] yields 
In 2). From this simple observation we infer that the approximation q <^ E used in 
the BT calculation is legitimate only at logarithmic accuracy, but not to calculate 
the constant term in ( p.. 21 ). The evaluation of this constant requires an accurate 
treatment of the very hard region q ~ Q'max- For the t-channel contribution we will 



derive the following analytic result correcting ( |1.2D , 



dx 



487r 



ln^ + ln24-7+^'^^^ 



C(2) 



487r 



In + 1.281 



2.031 



;i.3) 



where 7 ~ 0.577 is Euler's constant. 

We stressed that working beyond logarithmic accuracy requires considering the 
kinematic region of very hard transfers q E. Since k' = k + q and k < O (T), this 
corresponds to k' ~ E. Using 



s = {P' + K'y = 2P' -K' + ^ 2E'k'{l - cos^p/, 



;i-4) 



we infer from s ~ ET and k' E that the angle 9pik' between p' and k' must be 
small, and that the constant next to the leading logarithm (partly) arises from the 
angular region 

(1.5) 



Here is an essential point. When E ^ 00, the constant is sensitive to scatterings 
where p' and k' are coUinear. But collinear outgoing particles should be associated 
within the same 'jet' of particles, and such collinear configurations should thus be 
removed from the definition of observable collisional energy loss. Hence, we conclude 
that when E 00, the constant next to the leading logarithm in —dE/dx is not 
an observable, only the leading logarithmic term of ( |1.3D is meaningful. Strictly 
speaking, in the asymptotic limit E ^ 00 we can only state 



dE, 



dx 



4877 



In 



E 



0{l) 



1.6) 



where the constant ~ O (1) depends on the details of the jet definition. 

For finite E, however, 9pik' ~ \/T jE is a non-zero angle^. Depending on the ex- 
perimental angular resolution, processes transferring a large fraction of the incoming 
energy E to the particle K' (forming with P' the angle Q^i-ki 7^ 0) might be counted 
as observable energy loss. In the present study we assume the angular resolution to 



*For instance, for E = \Q GeV and T = 500 MeV we have ^TjE ~ 0.22 rad ~ 13° 
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be much better than ^JtJE. Then it is meaningful to include the (correct) constant 
next to the leading logarithm, as done in ( p..3|) . 

Comparing now ( p. .21) and (|1.3|) , it might seem that the difference in the constant 



term is only of minor importance. However, the incorrect approximation of Ref. [15] 
has also been used to calculate the coUisional energy loss of a heavy quark in QCD 
[9], and the latter must thus also be corrected. Moreover, the relevance of the 
region g ~ E' in the t-channel contribution suggests that s and -u-channels (Compton 
scattering) are important, contrary to what is assumed in [15]. Indeed, we find in 
the present case of QED (see section 3) that the Compton scattering contribution 
precisely arises from the domain ( |1.5|) ^. Since quasi-collinear configurations (|1.5| ) 
already contribute to the t-channel contribution ( |1.3D , Compton scattering cannot 
be dropped by invoking some collinearity argument. In fact, our calculation reveals 
that Compton scattering not only contributes to a new constant, but more crucially 
to an additional logarithmic term. 



dx 



967r 



In 



ATE 



M2 6 ^ C(2) 



;i-7) 



This logarithm is of collinear origin^, as is obvious from its divergence in the formal 
M ^ limit. 

In summary, if we aim to control the constant next to the leading logarithm in 
the t-channel contribution ( |1.3|) , s and u-channels must be included for consistency. 
This, in turn, brings an additional potentially large logarithm. For definiteness we 
state our complete result for the muon energy loss in the v —>■ 1 limit, obtained by 
adding to ([T3), 



dE 
dx 



4871 
487r 



2E ITE ^ 3 /C 2 

n — + - n — + n48 + - 
e^T 2 M2 2 V C(2) 

QfE 1 TE 
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This corrects the BT result ( p..2[ ). 

Our paper is organized as follows. In section 2 we focus on the contribution 
—dEf,/dx from t-channel exchange to the muon energy loss. In section 2.1 we present 
an exact relation between —dEe/dx and the muon self-energy. In section 2.2 we 
argue that the phase space can be conveniently decomposed in terms of the Lorentz 
invariant momentum exchange t, into the regions \t\ < \t*\ and > \t*\, where the 
cut-off t* satisfies e^T^ <C |t*| <^ T^, but otherwise is arbitrary. The contributions 



^The presence of s and u-channel contributions arising from the region q ^ E will also affect the 
QCD results obtained in [9]. 

^Such potentially large logarithms were mentioned in Ref. [18] in the case of the coUisional energy 
loss of light partons. 
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from the two regions are evaluated in sections 2.3 and 2.4 and summed in section 
2.5, where we derive the result (|1.3| ) and discuss the incorrect assumption made 
in the BT calculation [15]. (For completeness, we repeat the calculation of ( |1.3| ) in 
Appendices B and C by following the procedure used in Ref. [15], i.e. by decomposing 
the phase space into a soft and a hard domain with respect to q* = We can 

thus precisely see where the approximation q E used in [15] fails when evaluating 
the hard contribution from q > q*.) Section 3 is devoted to the contribution from 
Compton scattering —dE^/dx quoted in (|1.7|) , and our results are summarized in 
section 4. 

2. Scattering off electrons 

2.1 An exact relation: energy loss from self-energy 

Among the two processes contributing to the muon energy loss, namely scattering 
off thermal photons and scattering off thermal electrons (or positrons), we focus here 
on the latter mechanism, even though the idea of calculating the energy loss of a 
test particle from its self-energy is more general. Our discussion below follows the 
BT calculation [15] of the soft contribution to the energy loss, which we generalize 
appropriately. 

Let us start by recalling that the coUisional energy loss of a test particle is closely 
related to its interaction rate F, the latter being obtained from the imaginary part 
of the particle's self-energy E evaluated at the energy po = E + ie [19]. In the case 
of a muon, 

TiE) = (1 - npiE)) tr [(f + M) Im S(P)] . (2.1) 

For the t-channel interaction with electrons and positrons (involving a single photon 
exchange), the corresponding self-energy is the 1-loop graph depicted in Fig. |l|b, 
with a resummed photon propagator for reasons to be explained shortly. Since the 
on-shell self-energy is gauge invariant, one may choose a convenient gauge in order 
to evaluate (|2.1|). In Coulomb gauge, and without further approximation, the trace 




(a) (b) 

Figure 1: (a) t-channel scattering amplitude off electrons contributing to the muon inter- 
action rate, (b) The dressed muon self-energy. 
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m 



T]) reads [15] 

tr [(f + M) ImS(P)] = -47re^ (l + e"^/^) 



(icu (1 + nsiuj)) 



AB 

'2E' 



[2.2) 



A = Pl{uj, q) {2E^ -Eu- pq) + 2pt{uj, q) {p^ - Eu + pq - (pq/qf) 
B = {l-nF{E'))6{E-E' -uj)-nF{E')6{E + E' -uj), 



where E' = ^{p - qf + M^, Q = (cu, q) is the photon momentum, and Ub^f denote 
the Bose-Einstein and Fermi-Dirac thermal distributions. We also use the shorthand 
notation 

(2.3) 



(27r)3 

and the spectral functions of the longitudinal and transverse photons [15] 

TT 



(2.4) 



where /S.f t are the longitudinal and transverse photon propagators. 

According to the assumption M ^ T, we have npiE') <^ 1 and the factor B in 
^ reduces to 6{E — E' — uj). Performing the angular integral in ( |2.2| ) yields 



2ttv 



dq q 

+ pTi^^,q) 



duj{l + nB{uj)) <{ PL{i^,q) 

t UJ 



1-1^ 



t 



— + 



E 4^2 



Eq^ 2^2 {2Eqy 



■ (2.5) 



Comparing to the analogous BT result, we observe that relaxing their assumption 
|ti;|, g < T yields the additional terms put in between parentheses, and requires using 
the exact expression uj±{q) = E — a/ (p =F qY + M'^ for the bounds instead of the 
approximation ±vq. 

The energy loss per unit length is then obtained by weighting the differential 
interaction rate by uj/v, 



dE, 
dx 



2TtV' 



dq q 



dujuj{l + nB{uj)) \ PLi^,q) 



1-1^ 



t 



E 4^2 



+ Pt{uj, q) 



9 UJ 



tuj 



Eq^ 2^2 {2Eqy 



[2.6) 



As discussed in Ref. [15], the additional factor of u is improving the infrared behavior 
of the integral (|2.6| ) compared to (|2.5|) . In fact, evaluating Fg from Eq. (p.5|) with 
spectral functions obtained as the discontinuity of the 1-loop resummed propagators 
^i-ioop _ _ j|i-ioop^~i^ yields an infrared divergence due to soft transverse 

exchanges. In contrast, weighting the integrand of Fg by u leads to a well-defined 
energy loss 
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2.2 Phase space decomposition 

The necessity of expressing the energy loss in terms of a dressed exchanged pho- 
ton propagator (obtained by resumming the 1-Ioop photon self-energy) arises from 
the long-range nature of the Coulomb interaction. However, this infrared dynam- 
ics becomes unimportant for short-range interactions. This motivates, on physical 
grounds, a decomposition of the phase space in Eq. ( p.6|) . 

In the BT calculation of the energy loss, the matching procedure developed by 
Braaten and Yuan [20] is used. A momentum scale q*, chosen as eT ^ ^ T but 
otherwise arbitrary, is introduced to separate soft interactions, with q < q*, from hard 
ones with q > q*, see Fig. The soft contribution to dE/dx is evaluated using the 




q Q \f\ Q 

(a) (b) 

Figure 2: Two ways of decomposing the exchanged photon phase space, (a) Braaten- 
Yuan prescription, by introducing a cut-off q* separating a soft region (shaded) from a 
hard one, as applied in [15]. (b) Our approach, using a cut-off t* = — (f/*)^ in the invariant 
momentum transfer. 

'hard thermal loop' (HTL) approximation [21,22] for the dressed photon propagator. 
On the other hand, the hard q > q* contribution in (|2.6| ) is obtained by keeping only 
the leading term in the expansion p oc Im [(Aq ^ — H)"^] = Aq Im [11] -|- 0{ll^). This 
approximation corresponds to evaluating the elastic scattering amplitude in Fig. |l|a 
with the tree-level photon propagator (which would yield an infrared divergent result 
for the energy loss in the absence of the cut-off q*). The sum of the soft and hard 
contributions should, of course, be independent of the arbitrary scale q*, as was 
verified in [15]. However, this consistency check could not reveal that BT's result for 
the hard contribution is incomplete beyond logarithmic accuracy, as we will show in 
the following. 

For this purpose it will be convenient to decompose the phase space with respect 
to a cut-off t* = —{q*y in the invariant momentum transfer, as illustrated in Fig. ^3. 
This choice is motivated by two facts. First, in the region where \t\ < \t*\, the HTL 
approximation for the photon propagator is known to be valid [23] although u and 
q can be individually large. In fact, while the HTL approximation is usually derived 
under the assumption uj,q <^T (implying uj"^ + q^ <^ T^), it actually holds if the 
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Minkowski 'norm' \uj'^ ~ is small compared to [23]. Secondly, the calculation 
of the contribution where |t| > is more transparent (see section |2l^ ) , since in this 
region the squared scattering amplitude is a function of the Mandelstam invariants 
only. We stress that when |t| > \t*\ we can indeed neglect the 'medium modifications' 
to the matrix elements since the region \t\ > \t*\ is contained in the hard q > q* 
Braaten-Yuan region. 



2.3 Contribution from \t\ < 

We calculate this contribution to dEf./dx from (|2.6|) by changing variables to t 



q^ and x = co/q. Since \t\ < we can omit the terms of order Odt]^^"^ / E) 



and approximate the bounds on x by ±f . The term ^ uj/E in the integrand of ( |2.6| ) 
is easily checked to be exponentially suppressed. Using the HTL approximation for 
the spectral function we obtain 



dE, 



dx 



\t\<\t* 



dt {-t) 



dx 



X 



(1 -X2)2 



;i + nB(cu)) [pl + {v^ -x^)pt\ 



(2.7) 

where u = Xa/— — x^). From a simple parity argument, we can replace the 
factor 1 + nsiuj) by its even part, i.e. 1 + nsioj) |, and we find 



dE, 



dx 



\t\<\t*\ 



dx ■ 



X 



dt i-t) [pL + {v^ - x'')pt] 



Using the HTL longitudinal and transverse photon propagators 

1 ... 1 



(2.8) 



AL{uj,q) = 
with the self-energies [22] 

IlLix) 



q^ + Hiix) 



AT{uJ,q) 



^2 _ q2 _ nr(x) 



:2.9) 



ml [1 



Q{x)] 



^Ax) = [x^ + (1 - x^)Qix)] =^xil- x^)Q'ix) 



(2.10) 



Q(x) = — In , 

where m£, = eT/v^ is the Debye mass in the QED plasma, the expression 
becomes 



dE, 



dx 



\t\<\t-\ 



e2 


nv 


dx 


X 


87r2t;2 ^ 


1 ~v 


1 — x2 


e2 


I'V 


dx 


X 


87r2i;2 ^ 


1 —V 


1 — x2 



Im 



dt Im 



Hi In 



t 



t'2 — X^ 



t 



t-Urix) 



\t*\ + Ul - x^ 



1 — x2 t — nT'(x^ 

t1 + nT 



n, 



x^ 



Ur In 



(2.11) 
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We used ( p^) and introduced 11^(0;) = (1 — x^)IIl{x). We now take advantage of 
\t*\ ^ \IIl{x)\, \IIt{x)\, and extract the |t*|-dependence, writing for instance 



In 



\t*\ 



In- 



2\t* 



+ In 



m 



D 



m 



D 



We arrive at 



dx 



\t\<\t*\ 



487r2w2 



1 - 



2v 



In 



dx ■ 



X 



l + v 
l-v 

-Im 



In 



(2.12) 



2W 



m 



D 



x^ 



X 



x^ 



It In h 



(2.13) 



where 



(2.14) 



/i(x) = 2 ni(x)/m^ = 2(1 - x') (1 - Q{x)) 
frix) = 2 nr(x)/m^ = x(l - x'^)Q'{x) . 

The integral in ( p.l3|) seems difficult to evaluate analytically for arbitrary v. However, 
in the limit f — > 1 we are interested in, the result is surprisingly simple: the integral 
vanishes, as can be checked numerically. Hence, 

dE, 



dx 



m 



\t\<\t* 



(2.15) 



We end this section by presenting an alternative and fully analytical way to obtain 
the result (|2.15| ) - which proves indirectly that the integral in (|2.13| ) indeed vanishes 
for V = 1. Referring to Fig. ^, the kinematic region |t| < \t*\ is obviously given 



by the reunion of the soft region < {q 



*\2 



\t*\ and the region where q > \t* 



and \t\ < \t*\, or equivalently (1 — x^)|t*| < \t\ < \t*\, since \t\ = (1 — x^)g^. The 
contribution from the latter region is easy to calculate along the lines which led to 
( p.ll|) . When f = 1 we get 



dE, 



dx 



\t\<\t*\,q>q* 



87r2 
487r 



dx ■ 



X 



1 



x^ 



1 

1 3a;2 
dx In 

-1 2 



dt Im 



(l-x2)t* 
1 

— 



H, 



t-U 

487r 



t-Hr 
- In4 



(2.16) 



where we again used ^ |Hi(a;)|, |HT(a;)|. The contribution from q^ < \t* 
precisely the BT 'soft' contribution, which reads for f = 1 [15] 



IS 



dE, 



dx 



q<q* 



247r 



In + 0.256 

eT 



(2.17) 



In Appendix B we rederive this result and determine the constant analytically. 



dE, 



dx 



q<q* 



247r 



, q* In 24 
n — H 

eT 2 



(2.18) 



By adding ( |2.16| ) and (|2.18| ) we confirm the result (|2.15|) for the contribution from 
the domain \t\ < \t*\. 
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2.4 Contribution from \t\ > \t* 



As already argued in section |2.2| , wlien |t| > |t*| thermal corrections to the exchanged 
photon propagator can be ignored. In this kinematic domain the energy loss can thus 
be obtained from the general relation [15] 

dEi _ 1 friiik) f ni{k') f 1 ,^„,4x(4)^D , rWNl>v-iA^i2^ 



dx 2Ev 



A 2k A, 2k' J^,2E'^ ' ^ U^J 



spins 

(2.19) 



by inserting a factor 0(|t| — |t*|) in the integrand. The tree-level matrix element 
TVlj describes the scattering off a target particle of type i. Summing over 
initial and final spin states and dividing by the degeneracy factor d of the incoming 
test particle gives its energy loss —dEi/dx averaged over spin states. Furthermore, 
ni{k) = {exp{k/T) ± 1)^^ is the thermal distribution of the target particles, and 
•hi = 1 ± Hi accounts for the Bose enhancement or Pauli blocking for the scattered 
state. In line with our previous considerations, the target particles are assumed to 
be massless. 

The tree-level matrix elements squared depend only on the Mandelstam inva- 
riants s and t. Then, as derived in Appendix A, the phase space integral in ( p.l9| ) 
can be reduced to 

dEi 1 , f nAk) f s + k\ r° , , , dai 

- dJ^ l-——--] dti-t)-^, (2.20) 



dx V J. 2k \ s - M"^ E J Jf . ^ ' dt 

where only the approximation 1 has been made, which is also justified in Ap- 

pendix A. Apart from this simplification, the expression (p.20| ) is exact, in particular 



with respect to the kinematics of the scattering process. The bound in the t-integral 
is 

_ {s-M^f 

^min ; \2-2^) 

S 



and we introduced the differential cross section 

dai 1 



dt lQ7r{s- M^y dd 



' (2-22) 



spins 



where di is the spin degeneracy of the target particles. 

Our expression ( |2.20|) generalizes a formula used by Bjorken [1] for a massless 



test particle of infinite energy, to the massive and finite-energy case. It will allow us 
to calculate in a rather compact way the contribution from |t| > \t*\ to the coUisional 
energy loss. Focusing on the large E limit and inserting 0(|t| — |t*|) in the integrand 
of (|2.20| ) we obtain 



dEi 



dx 



\t\>\t*\ 
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The energy loss —dEg/dx arising from \t\ > \t*\ of a muon scattering off electrons 
and positrons'' is obtained from (|2.23| ) with (|2.22|) and 



2 5^|M 



spins 



[s-MY s i 
— H h - 

t2 ^ t 2 



(2.24) 



We mention that the two last terms of ( p.24|) contribute to the t-integral in ( p.23| ) as 
a constant, 



1 



t* 



dt (-t) 



s 1 
1+2 



3 

'4 ' 



(2.25) 



where we used s 3> ^ \t*\. This results in a non-zero contribution to —dE^/dx, 



487r 



(2.26) 



which arises from the very hard region —t ~ — tmin ~ i.e. q ^ O [E). This proves 
that the assumption q <^ E used in [15] is inappropriate when calculating —dE/dx 
beyond logarithmic accuracy. 

Only the first term of ( |2.24|) is sensitive to the cut-off t*, since for t* — > it yields 
a logarithmic divergence in the infrared, 



1 







f dt {-t) 


t^ 



With the help of the definite integrals 



In — — — ~ 



In 



(2.27) 



dx npix) X 
dx npix) xlnx 



71 

12 
12 



1 - 7 + ln2 + 



m 

C(2), 



(2.28) 
(2.29) 



the contribution of the first term of ( ^.24| ) to ( ^.23] ) is found to be 



8TE C'(2) 
In -r— ^ 7 + 



t* 



C(2) 



(2.30) 



Adding the contributions ( p.26|) and (|2.30|) we obtain 



dE, 



dx 



\t\>\t*\ 



487r 



■ SET 3 C'(2) 

In -; — : 7 ' 



t* 



4 C(2) 



(2.31) 



''The contribution of positrons is identical to that of electrons for an e*7 plasma with vanishing 



chemical potential, and is accounted for by a factor 2 in Eq. (2.24) 
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2.5 Complete t-channel result and discussion 



Combining the contributions from |t| < \t*\ and \t\ > \t*\ given by ( |2.15D and 
we find the muon energy loss from scattering off electrons and positrons as 



dx 



487r 



ln^ + ln24-7+^'^^^ 



C(2) 



(2.32) 



as already quoted in ( |L3|) . 

We stress that the leading logarithmic term in 



( ^■32| ) arises from the region 



< -t < s ~ ET. Since ~t = -{K - K'Y = 2kk'{l - cos Okk') and A: ~ T, this 
implies k' = \k + q\ <^ E. The approximation q <^ E used in [15] is thus legitimate at 
logarithmic accuracy. Beyond logarithmic accuracy however, the region contributing 
to ( p.32| ) extends to mj^ ^ ~t ^ s ^ ET. In particular the region of maximal 
transfers —t ~ — tmin ~ s ~ ET {i.e., k' ~ E) affects the constant. Thus the latter 
could not be correctly determined in [15], where the approximation q <^ E was used. 
As mentioned in the Introduction (see (|1.4|) and ( |1.5| )), k' E also implies that the 
angular domain Opik' ~ \JT/E -C 1 contributes to the constant. 

The constant in our result ( p.32|) differs from that of the BT result ( |1.2|) . In order 
to confirm our result, we present an alternative calculation in Appendices B and C. 
There we follow the BT approach by using the familiar Braaten-Yuan decomposition 
of the phase space into a soft q < q* and a hard q> q* kinematic domain, see Fig. 
In Appendix B we confirm (see ( |B.12|) ) the BT result ( 2.17 ) for the soft contribution 
in the f — ^ 1 limit. In Appendix C, we repeat the BT calculation of the hard q > q* 
contribution without using the approximation q <^ E. As expected, we find that the 
constant is sensitive to the hard domain q E. As shown explicitly in Appendix C, 
the approximation q <^ E would amount to neglect terms both in |A^e|^ and in the 
5-function for energy conservation, which are important to determine the constant. 
Those corrections lead to the result ( |C.24| ), instead of ( |C.23| ) as found by BT^. 

Adding the soft ( p.l2|) and hard ( p.24|) contributions we recover ( p.32|) , found 
within our decomposition of phase space using an invariant separation scale t* . We 
view this as a corroborating evidence of the correctness of our results. 



3. Compton scattering 



The contribution —dE^/dx to the muon energy loss from Compton scattering (see 
Fig. is dominated, as we will shortly see, by hard transfers —t ~ s. Thus it can 



^It is a coincidence that the last term ~ —3/4 of ( C.24 ) missed by BT is identical to the 



contribution (2.26). Hence, in effect the BT result is the same as what we would obtain by keeping 
only the term - (s - APf /t^ in (p. 241). 
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be obtained from ( p.23|) by setting t* = and by using 



8e^ 



spins 



-U 



2M- 



— + 2M2 

-u 



+ 2M- 



2 5 + u + 
su 



(3.1) 

where we define u = u — and 5 = s — M^, which satisfy s + u + t = 0. The 
three terms in (|3.1|) correspond to the contributions from the s and w-channels, and 
from the interference term. Using (|3.1| ) in our formula ( 2.23| ) we can easily show that 





Figure 3: Amplitude for Compton scattering. 



all terms which explicitly depend on in ( |3.1| ) yield contributions to —dE^/dx 
which are suppressed by at least one power of s ~ ET when E ^ oo. Thus for our 
purposes ( |3.1| ) can be approximated by 



8e' 



spins 



— + 

S 



-u 



This yields the integral 



57r 



,^u + s 
du 



—u s 
s —u 



+ 0{l/s) 



In 



M2 6 



+ 0{l/s) 



(3.2) 

(3.3) 
(3.4) 



where we changed variables from t to -u, with the bounds -Umin = —5, ttmax = — s/ s. 

The leading logarithm arises from the w-channel, more specifically from the kine- 
matic region u^i^ ^ tt ^ "Umax, i-G- M"^ ^ ~u ^ s when s ^ M^. Since in this 
region — t ~ s ~ s (recall that s + u + t = 0), the physical interpretation of the 
logarithmic enhancement is the same as for the total cross section cr^ = j dt da^/ dt. 
The latter behaves in the high energy limit as a"^ ln{s / M"^) [24], with the loga- 
rithm originating from backward scattering of the photons in the center of momentum 
frame. We thus infer that the Compton contribution to the energy loss is dominated, 
at large energies and to leading logarithmic accuracy, by the same mechanism. 

Plugging (|3.4|) into ( p. 231) , and using 

^2 



77 



dxnsix) X = — 



dxuBix) xlnx 



TT 

~6 



1-7 + 



C(2)J 



(3.5) 
(3.6) 
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for the integral over k, we obtain 



dx 



967r 



In 



ATE 5 C'(2) 

_L 2_L^ 

M2 6 ^ C(2) 



(3.7) 



as quoted in (p..7|). 

The logarithm in ( ^.T] ) arises from ^ — tt ^ s, while the regions —u ~ 
and —u ~ s only contribute to the constant. Hence the complete expression ( |3.7| ) 
stems from^ < — n < s, including also very hard exchanges — t ~ s ~ ET ^ M^. 
/,From — t = —(-ft' — K'Y = 2kk'{l — cosOkk') ~ s ~ ET, the typical values of 
k' = k + u! contributing to (|3.7|) are k' ~ ^max ~ ^- From (|1.4|) we thus find that s 
and ti-channel contributions arise from the angular domain 6p'k' ~ a/TZ-E ^ 1, as 
anticipated in the Introduction. 

It is instructive to write the integral ~ / du/u appearing in ( p^) as an integral 
over the angle between p and k' by using 



u 



{P - K'Y -M^ = -2PK' ~ -2Ek'{l - cos^ 



pk' 



ET, we must have 6pk> ^ 1 in 



(3. 



and the 



Hence, in order to satisfy —u < 
domain < —u < s leading to the logarithm translates into M/ E < Op^i < ^JT jE. 



(3.9) 



Consequently, the result ( p.TD stems from the angular regions 

MIE < V < Vt/e, ep>,> ~ Vt/e. 



4. Summary and outlook 

In this study we have reconsidered the Braaten-Thoma calculation [15] of the muon 
collisional energy loss in a hot QED plasma. For the t-channel contribution (scatter- 
ing off electrons and positrons), we have shown that controlling the constant next to 
the leading logarithm requires, in the limit E ^ M'^/T, considering the region where 
the invariant transfer —t is on the order of its maximal value, —t ~ — tmin — s ^ M^. 
The BT calculation of the t-channel contribution, which was based on the invalid 
assumption q <^ E, evaluated the constant incorrectly. We obtained the corrected 
result for this contribution in Eq. ( |1.3| ). 

We showed that the 'constant' is sensitive to the angular domain Opiy ~ ^JtJE, 
i.e. to coUinear configurations when E —>■ oo. Thus, for consistency the contribution 
from s and w-channels (Compton scattering), which arises from similar configura- 
tions, must be included in the energy loss. As already stressed in the Introduction, 
this is our main message. Removing Compton scattering from the definition of en- 
ergy loss implies that we have to give up determining the constant next to the leading 

^Since M ^ T, the contribution from s and u-channels arises from exchanges where thermal 
corrections to the muon propagator are suppressed. The calculation using a bare muon propagator 
in Fig. is thus legitimate. 
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logarithm in the t-channel contribution. In other words, working beyond logarith- 
mic accuracy is meaningful only with Compton scattering taken into account. The 
Compton process yields a potentially large 'coUinear' logarithm oc InlET/M"^), see 
Eq. (|1 . 7|) . It arises from hard transfers —t ~ s, and was previously neglected in [15]. 

It will be interesting to study the consequences of our findings for the collisional 
energy loss of a heavy quark in a hot QCD plasma, and to see how the results of [9] 
are modified^". We also point to the necessity, in phenomenological studies, to take 
into account the finite experimental angular resolution 6. If the latter is of the order 
of \/T/E (or larger), the final state configurations with 6p'k' < S should be removed 
from the definition of energy loss, modifying our full result (|1.8| ) by introducing a 
(5-dependence. Since (quasi)-collinear configurations 9p'k' ~ \^^TJE correspond to 
hard exchanges q E, the angular resolution 6 will actually translate to an upper 
cut-off in q, above which the elastic processes under consideration will not contribute 
to an observable energy loss. 

A. Thermal phase space 

For a given function f{s,t,uj) depending on the Mandelstam invariants and the 
energy transfer u = E — E' , we calculate the functional 

I[f] = ± f^Kk^ f n{k;i f J_^2n)'s('\P + K-P'-K')fis,t,u;). (A.l) 
^ 2E 2k J I,, 2k' Jp,2E'^ ' ^ ) J\ ^ ^ > \ ) 

We start with the /c'-integral, for which the specific form of f{s,t,uj) is not 
relevant, since s and t are determined by k, p' and p only. Following a standard 
procedure we write 

and evaluate 

i ^ (2^)'^^'^(^ + K-P'-K')= 2nn{K)e{K) S (^'^) . (A.2) 

Here = K+P—P' is fixed by momentum conservation. With P—P' = Q = {uj,q), 
we have in particular /cq = k + u. 

In order to proceed with the p'-integral, we specify a coordinate system. We 
choose the 2;-axis along the direction of p, and orient the yz-pla.ne to contain k, 

p= (0,0,l)p, 

k = {O,sinip,cosip)k , 

p' = (sin6' sin 0, sin6' COS0, cos6')p' . (A. 3) 

^"A first step in this direction was done in [25]. 
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The integral over the azimuthal angle (j) is readily performed with the help of the 
^-function in (|A.2| ). We first express its argument in terms of the Mandelstam in- 
variants t = and s = (ir + Pf = + 2KP, 

^'2 = {K + Qf = 2KQ + Q^ = 2K{P - P') + t = s - M"^ + t - 2KP' . (A.4) 

Writing KP' = kE' — kp' and using ( |A.3| ), we find K '^ = A + B cos (p, with 

A = s- M"^ + t- 2kE' + 2kp' cos ipcosO, 

B = 2kp' sin ipsinO . (A. 5) 

Consequently, 

/■27r r, 

/ d<p6{K:')=—e{g), (A.6) 
Jo V9 

with g = B^ — A^. The B-function reflects the kinematic constraints imposed by 

energy-momentum conservation. 

For the remaining integrals we change variables, from p' and cos 6 to 

t = 2 {M^ - EE' + pp' cose) , 

uj = E-E', (A.7) 
with the Jacobian E'/{2pp''^). Using ( |A.2|) , ( [A.6|) , the expression ( |A.lj ) becomes 

As mentioned above the precise kinematic bounds on the u and t integrals will 
naturally arise from the condition g = B^ — > 0. We easily obtain from ( [A .51) 

g{uj) = -a'^uj^ + buj + c, (A. 9) 

whose coefficients can be expressed as 
s- 

a = , 

p 

9f 

6= {E{s-M^)-k{s + M^)), 

c=-^[t {{E + kf -s)+ Ap^e -{s-M^- 2Eky] . (A.IO) 

Because the quadratic term in ( [A .91) is manifestly negative, g{uj) is positive in an 
interval [uJmin, t^max] where the discriminant D = Aa?c + IP' is positive. We have 

, .max _ b±VD 

^min - ^ ) {-^■^^) 

^ / / 9n9\ /4fcsin'?/'\^ 

D = -t{st + {s-M^f) i — ——] . (A.12) 
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The condition D > leads to the famihar range of the invariant momentum exchange 
in 2 ^ 2 processes with one massless and one massive colhsion partner, namely 
^min ^ t < with 

^min = • (A. 13) 

S 

We now show that when g{uj) > (or equivalently cUmin < < t^max), the factor 
Q{k+u!) appearing in ( |A.8| ) is actually redundant. Recall that the condition g^u) > 



arises from ( A. 6 ). Thus the values of u contributing to the l.h.s. of (|A.6|) , i.e. for 
which K '^ = 2KQ + t = 0, must belong to the interval [comin, i^max]- Such values thus 
satisfy 

t + 2kuj = 2kq {t + 2kujf < Ak'^q^ = 4k^{uj^-t) t{t + 4k{k + uj)) < . (A.14) 



Since t < we find that k + u > -t/{4:k) > 0. Thus the e(A; + u) factor in (^) 
can be dropped and we can now specify the precise bounds on t and uj: 

m = f dt n^^n{k + uj)f{sAu^). (A.15) 

As a side remark, let us note that the change of variables ( |A.7| ) maps the original 
integration area, p' G [0, oo[ and cos^^ G [—1, +1], into a (t,uj) region enclosed by 



^From this expression we check that the maximal energy transfer is as expected 
Max(u;+) = E — M, occurring at t = —2M{E — M) and corresponding to 'full 
stopping'. Note also that our derivation of the bounds on u in ( |A.15| ) implies that 

UJ-{t) < LOmin < UJ < UJ^s,^ < UJ+{t). 

Approximation n ^ 1 



The integrand in formula ( |A.15| ) contains the thermal distribution n{k+uj), which 



usually prevents the calculation of the cu-integral in terms of elementary functions. 
We may, however, obtain useful approximations of the integral by replacing n = 
1 ± n — i> 1, ie., by neglecting thermal effects on the final states. 

Under this assumption we can evaluate analytically the a;-integral in ( |A.15| ) for 
the function f{s,t,uj) = uj^ \A4\'^{s,t), for instance by using the formal identity 



UJ „ I , UJ 



/(f) = / duj^ = Re duj^ (A.17) 

Up to the prefactor Q{D) involving the discriminant D of the quadratic function 
g{uj), which reflects the 2-body kinematics and ensures a non-zero support of the 
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integral as discussed above, we obtain for example 








1 


2 




TT 

a 


71 b 
2 ^ 


7T D + 2b^ 
8 a5 



(A.18) 



Thus, replacing n by unity in ( |A.15|) , we find 

1 



I[f] ^ JW[A^] 



^1 (A.19) 



167r2pE 

For £ = 1, as of interest for the energy loss calculation, the approximation 
n = 1 ±n ^ 1 can be justified, as was done in Ref. [15], and as explained here below 
in Appendix C (after Eq. ( p.6|) ). We will use this approximation, which allows to 
write, using HKJEj ) for £ = 1 and (^jKJOj ), 



n{k) 
~2k~ 



s + M"^ k 
s-M'^E 



dt {-t) 



\M\\s,t) 
167r(s-M2)2 



(A.20) 



which we recognize as an integral of the differential cross section weighted by the 
factor t. 



B. Scattering off electrons: soft q < contribution 



We start from Eq. (40) of Ref. [15], which reads in the limit f — 1: 



dx 



soft 



i7T 



/ dq^ duJUJ [pL{uJ,q) + [1 - x'^)pT{uj,q)] , (B.l) 

Jo J-q 



where x = uj/q. For v = 1, the u integration range is the whole space- like region 
|a;| < q. We can shift to the time-like \uj\ > q region by writing 



duj 



doj 



doj . 



(B.2) 



u)\>q 



In ( B.l ) the resulting integral over the infinite cu-range vanishes, which can easily be 
seen from the sum-rules^^ 



dujuj Pt{uj, q) = I , 

/oo 
duu^ Pt{uj, q) = q^ + 
-oo 



oo 
oo 



m 



D 



duj uj pL^i^, q) 



m 



D 



3q^ 



(B.3) 



11 



The sum rules can be derived from the spectral representations of the gluon propagators [22] 



-At(u;, k) = I dk,P^^^ , -AUu^, k) 



1 

fc2 



dki 



PLiko,k) 



ko — Lu 



ko — Lu 

by identifying in the w ^ oo expansion of the left and right-hand sides the terms of appropriate 
order in 1/uj. 
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where m|, = 

In the time-hke region |ti;| > g, the spectral functions are given by the pole 
contributions (s = L, T) 

Psi^,(l)\\uj\>g = 4^)zs{q)S{uj^ -uj^,{q)) , (B.4) 

and ( |B.1|) becomes 



dx 



soft 



e 

'8^ 



dq^ [zL{q) + (1 - xl)zT{q)\ 



(B.5) 



where we denote Xs = ujs{q)/q- 

By definition, the poles = uJs{q) of the propagators ( |2.9| ) satisfy the implicit 
equations 

(B.6) 



q = -IIl{xl) , g - 



with the longitudinal and transverse photon self-energies as specified in ( p.lO|) . 
The residues Zs{q) are defined by 



Zs{q) 



Zs{q) 



uo"^ — ujg{q) q^{x'^ — xD 
Expanding the denominators in ( |2.9D around x"^ ~ x1 we get 



when ^ ^liq) = ^l{q)/q^ ■ (B-7) 



ZL{q) = 2xl — 
zriq) 



q 



U'^{xl) 

1 



-2a;, 



1 



1 - 



2xt TlTixT) 



(B.8) 



Thus zi^riq) are explicit functions of xl,t, which suggests to shift variables from 
to Xl,t in ( p.5|) . For the longitudinal and transverse contributions we find from ( p.6| ) 
and {JKW) 



ZL{q)dq^ = 2xLllL{xL)dxL 



[1 - XT)zT{q)dq^ 



2xt 



IlT{xT)dXT ' 



(B.9) 



^From (p.5|) we then obtain 



dE,, 



dx 



soft 



e 

4:71 



/ dxxUL{x)+ / dx ^ ' ; 

JxL{q*) JxT{q*) X —1 



(B.IO) 



Using now g* ^ eT we have XL^xiq*) — 1? and since the first term of ( p.lCI| ) is 
integrable at a; = 1, we can safely replace x^lq*) —>■ xxiq*) in this term. Using the 
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relation IIt{x) = n'^(x)x(a;^ — l)/2 (obtained from ( p.lOj )), the two terms of ( B.10| ) 
combine into a full derivative, to give 



dx 











soft 





D 



xriq*) 



4 1, 2 

h - In — 

3 2 XT{q*) - 1 



Using finally XT{q* ^ eT) ~ 1 + mjj/^Aq*^) [22] we arrive at 



dE, 



dx 



247r 



soft 



, In 24 4 
In — H 

eT 2 3 



(B.ll) 



(B.12) 



C. Scattering off electrons: hard q> contribution 

The hard contribution to {—dE/dx) reads [15] 



dE, 



dx 



_ I f I 

hard ^ J p' 
X 



np{k) f npik') 



E .L 2E' jf^ 2k y^, 2k' 



spins 



where uj = E — E' is the energy transferred by the muon in the elastic scattering. 

The squared t-channel scattering amplitude (summed and averaged over spins) 
is given by 



1 ^ = 16 ^ [{PK){P'K') + {PK'){P'K) - M^KK'] . (C.2) 



spins 



In BT it is assumed that k' ~ T, which allows for the approximation 



spins 



2{k - vk){k' - vk') + 



2^2 



(C.3) 



where v = p/E is the incoming muon velocity. However, ( C.2 ) can easily be cast in 
a form similar to (|C.3D without any approximation. Using P' = P + K — K' and 
KK' = —PQ = —t/2 we obtain from (|C.2|) the exact expression 



1 J2 |A<r = 16^[2{PK){PK') + {M^ + t/2)t/2] 



spins 



MH e 



2{k - vk){k' - vk') + -— + 



2^2 4^2 



(C.4) 



which differs from the BT approximation ( |C.3| ) by the term oc / E"^ and the pref actor 
(^2 instead of EE'). 
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Recalling that uj = k' — k, the (5-function for energy conservation in (C]T) can be 
expressed without approximation as 



6{E - E' -u) /{2E') =6{{E 



E'^) = 6 {2pq + t - 2Eu) 



6{u -vq-t/{2E)) /{2E) . 



(C.5) 



Using (|C.4|) and (|C.5|) , and 3-momentum conservation to perform the integral over 
p', the energy loss ( |U.1| ) becomes, again without any approximation 



dE, 



dx 



IQne 



hard 



2{k - vk)ik' - vk') + -— + 



2^2 4^2 



(C.6) 



This expression can be simplified as follows. First, as in the BT calculation, the 
term oc npik') is neglected, i.e. we replace fipik') — > 1. Indeed, for k' ^ T, npik') 
is exponentially suppressed. For k' ~ T, on the other hand, we have uj = k' — k T 
and g = — fc| ~ T, thus t/{2E) ~ T'^/E can be neglected in the 5-function for 
energy conservation, and the term oc np{k') then vanishes by antisymmetry in k ^ k' 
(recall that t = {K — K'Y). Secondly, in the integrand of ( |C6| ) we insert 



1 = j d^q5%q + k-k') J du;5{uj + k-k') 



(C.7) 



and perform the integral over k'. This gives 
dE, 



dx 



hard 



47re^ rnp{k) 
V Jk k 



duj 



5{uj + k- 


k + q\) 




k + 


q 





5{uj -vq-t/{2E)) 



t MH 



2{k - vky + {k- vk)— + + 



E 2^2 4^2 



(C. 



Eq. ( |C.8| ) differs from the BT expression only by the t/{2E) term in the 5-function 
for energy conservation and the t^/lAE"^) term in the expression of the squared am- 
plitude. 

As noted in BT, since —dE/dx does not depend on the direction of v, it is 
convenient to average ( |C.8| ) over this direction, using 



4^^^"-"^) = ' 

—5{uj - vq) = '—q' 

4tt 2vq q 



(C.9) 



— 6(lj - vq) f V = ^ — 

An ^ ^' 2vq 



— 5 + — qq^ 



2g2 



2g2 
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where we use the notation d) = lo — t/{2E). Using then 5{uj + k — \k + q\) 
2\k + q\5{t + 2kuj - 2kq) we find 



dx 



47re^ 



hard 



/ /- r'^'^ duj6{t + 2kuj-2kq) 

Jk ^ JqQ Jcu^iq) 



X 



kto v'^q^ — 



2\k' - 2kq— + 



■^[k-kq- 



g2 2g2 

-T - + — - + 



3cD2 - t;2g2 



gV ^ 2^2 4^2 



t2 



2g4 



(C.IO) 



where the factor 0(t> g — ) yields the bounds on tu, 



w±(g) = E - v/^2 + g2 ^ 2Evq . 



(C.ll) 



/^From ( |C.1CI| ) we proceed as follows. We replace 2kq ^ t + 2kuj in the integrand, 
perform the integral over the angle between k and q using 



dcoi^e5{t + 2ku-2kqcoi^d) = e{\q-k\<\u + k\<q + k)/{2kq) , (C.12) 



and then re-express a; — >• a; — t/{2E). Ordering in powers of we obtain 



dE, 



dx 



hard 



47r^f 2 



dknpik) I dq 



duQ{\q - k\ <\u + k\ <q + k) 



uj^(q) 



UJ \ ?>UJ 



g2 



4g2 



l-f^g^ kik + u) 



k{k + u) 



4 



t 



uj[12k{k + u) + 3uj^ -q^] , Ak{k + uj){3uj^ - q^) + 3uj^ -2u^(fj-q^ q^ \ 

^4E2J ■ 



Aq^E 



16^2^2 



(C.13) 



We mention here that all 1/E and l/E"^ terms stem from our 'correction' t/{2E) 
in the 5-function of ( p.8| ), except the last term ~ q'^/{4E'^), which comes from the 
exact expression ( |C.4| ) of the squared amplitude. We have separated these terms to 
underline (see ( |C21[ ) and ( |(J.22| )) that the error made in [15] is due to an incorrect 
approximation both in the 5-function for energy conservation and in the expression 
of the squared amplitude. 

Using —q < uj-{q) < uj+{q) < g we can show that the phase space constraints in 
( |C13| ) can be written as 



e(g - q*) 0(|g - k\ <\uj + k\ < q + k) e(cu_ <uj<uj+) 
Q{k - q") Q{q* <q<k) Q{u^ < lu < 
+e(Max(fc, q*)<q< gmax) Qiq - 2k < lu < iu+) , 



(C.14) 
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where k and gmax are the values of q at which q — 2k = uj-{q) and q — 2k = uj^{q), 

J2k{E + kl_ 2k{E + k) 

E{l+v) + 2k' ^'""^ ~ E{1 -v) + 2k- ^ 

At this point we use k ^ T <^ E to approximate k ~ 2k /{I + v), and q* <^ T so 
that effectively k > q*. Also, in the part of the integral oc G(g* < g < A;) we can 
approximate uj±{q) ^ ±vq. We can thus replace in ( |C.13| ) 



dq duj<d{ 



2k/{l+v) rvq /-gmax r^+{l) 

dq duj + dq duj 

* J-vq J2k/{l+v) Jq-2k 



(C.16) 



In the term corresponding to g* < g < 2k/{l + v), q and u are constrained to be of 
order T. It is then easy to see that in the curly bracket of ( p.l3|) , the first line will 
contribute as ~ e^T^ to the energy loss, whereas the second line can be neglected, 
since these terms are suppressed by 0{T/E) and 0{T'^/E'^), respectively. For the 
contribution from 2k/{l + v) < q < g^ax in ( |C16| ), the c<j-range brings a factor 
~ A; ~ T, as can be seen from the identity 

.,iq)-iq-2k) = —r^^^±^L__(i_JL \ . (C.17) 

' ^ ' 2{E + k)-q~uj+{q) \ q^^J ^ ' 

When g ~ g^ax ^ we have there c<j ~ g ~ gmax, and the second line of the 
curly bracket of ( p.l3| ) contributes to O {qmax/ E) and O {q"^^^/ E^) (the terms oc k 
are suppressed by at least O {T / E) and can be dropped). This contribution is thus 
important when g^ax ~ E. From ( |(J.15| ) this happens when £'(1 — v) < k ^ T, i.e. 
when E > M^/T. 

Using ( |C.16|) the expression (|C.13| ) can be written as 



dE, 



dx 



hard 47r3t;2 



dknpik) 



2k/{l+v) .vq r-^+ii) 

— duLU + — duouo 

q J~vq J2k/{l+v) q Jq-2k 



3uj^ 1-v^q^ kik + u) ,^ ^^k(k + uj) 



e 



4g2 4 ' 2 Q2 ' q2 ' ^ . 

4 /""^i. _fi,^ T"^'' dq r+^'^, f uj{q^-3uj^) 3uj^ - 2uj^q^ - q^ q^ \ 
2k/ii+.)7 Jq-2k M'E + We^ +4E^/- 

(C.18) 



This expression differs from the BT pendant by the upper bound uj^{q) (instead 
of vq) in the first line, and most importantly by the presence of the third line, 
containing terms formally (x 1/E and oc l/E"^. In order to compare further to the 
BT calculation, let us consider as in [15] the two limiting cases E ^ M'^/T and 
E ^ M'^/T, where the expression (|C.15| ) of gmax can be approximated as 

2k E"^ 

gmax ^ , Z < E VS. gmax - E. (C.19) 

E«Ml 1 - V s»m£ 
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In the domain E <^ M'^/T, we have q < Qn 
For E 3> M'^/T, we can approximate 



<^ E, thus uj+{q) ~ vq from ( |C.11|) . 



= E - + ^2 _ 2Et;g = E - W(E - g)2 + ^ 

V 1 + f -C/ 



(C.20) 



where we assumed E — q ^ M. This is justified since the contribution from E — 
M < q < E to the energy loss ( |C18| ) is of order e^T'^M/E and thus suppressed 
compared to the dominant contribution ~ e'^T^ we are looking for. Thus either 
when E -C M^/T, or when E ^ M'^/T (corresponding to the ultrarelativistic v —* 1 
limit), the approximation ~ vq in ( P.18|) is valid^^, and the only difference 

between our result and the BT calculation is the additional term written in the third 
line of ( |C.18|) . As discussed previously, this term contributes to the energy loss as 
~ e^T^ only when E > M^/T. We conclude that the domain E > M^/T {i.e. 
V 1) is treated incorrectly in [15]. Focusing now on this limit, we obtain from 

(EH) 



dE,, 



dx 



dE,, 



hard 



dx 



hard,BT 



47]-3 



oo /-E 

dknpik) 

Jk 



dujuj 



Q Jq~2k 



2E 



+ 



4^2 j ' 
(C.21) 



where we replaced u —>■ q in the bracket of the uj integral, the terms oc {q — u) 
yielding negligible contributions when E ^ oo. As mentioned after ( |U.13| ), the 
corrections to the BT result arise from using the exact squared amplitude (|C.4|) (term 
~ q'^/{4E'^) in ( |C.21| )) and the (5-function for exact energy conservation (additional 
term ~ —q/ {2E)). The remaining integrals in ( P.21|) are trivial and we get to leading 
order in 1/ii^ 



dE, 



dx 



dE,, 



hard 



dx 



+ 



hard,BT 



487r 



(C.22) 



Our new term in ( |C.22|) arises from a kinematical domain where the momentum 
exchange q is 'very hard', q ~ gmax — E. This domain already contributed to the BT 
result (written in ( |C23D below). Indeed, the logarithmic term oc In E/T arises from 
an integral ~ J'^ dq/q, where, for instance, the interval E/2 < q < E contributes as 
In 2. The very hard region was however not consistently treated in [15], due to the 
ad hoc use of the approximation q <^ E. 
Given that [15] 



dE, 



dx 



hard,BT 



487r 



In 



2TE 8 



7 + 



C(2)J 



(C.23) 



^This approximation would be incorrect in the intermediate regime E ~ IvP /T , where the exact 



expression (C.ll) of 0J+{q) should be used. 
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our result ( C.22 ) reads 



dx 
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